Gluon-Glueball Duality and Glueball Searches 
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We discuss a notion of gluon-glueball duality analogous to quark-hadron duality. We apply this 
idea to the radiative decay of heavy orthoquarkonium, QQ — > ■ygg, which has been used to search 
for glueballs. The duality is first introduced in two simplified contexts: (i) a hypothetical version 
of QCD without any light quarks and (ii) QCD in the large- N c limit. We then discuss how an 
approximate form of this duality could hold in real QCD, based on a hierarchy of time scales in 
the temporal evolution of the gg subsystem in radiative orthoquarkonium decay. We apply this 
notion of gluon-glueball duality to suggest a method that could be useful in experimental searches 
for glueballs. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) is very successful 
theory describing quark and gluon interactions. There 
are ample observations of gluon jets in high-energy col- 
lider data, and lattice QCD calculations of the pure glu- 
onic sector of the theory have yielded a detailed spectrum 
of (color-singlet) bound states of gluons, commonly called 
glueballs [IJ-Q- The lightest of these can be modelled as 
gg states, where g denotes a gluon; these include a state 
with S = 0, L = 0, and J PC = ++ , and a heavier state 
with S = 2, L = 0, and J PC = 2 ++ . Radial excitations, 
states with angular momentum L > 1, and ggg states 
also appear in the spectrum. Over the years there have 
been numerous experimental searches for glueballs. It 
was pointed out early on that a promising method is to 
use the radiative decay of a heavy QQ orthoquarkonium 
state [4, 5]. At the level of elementary constituents, this 
decay is QQ — > jgg, so that when the two gluons are 
emitted with an invariant mass close to that of a glue- 
ball, they have substantial probability to bind to form 
this state. Other production channels have also been 
used. At present, there are strong indications for hadrons 
with large gluonic components, although there are still is 
no consensus concerning the details of the mixing of qq 
and gluonic components to form various physical mass 
eigenstates 

In this paper we examine the temporal evolution of 
glueball production in radiative orthoquarkonium decay. 
We use the fact that glueballs have a smaller density of 
states than qq mesons, as a function of mass, in con- 
junction with the Heisenberg uncertainty principle, to 
infer that one can generically measure the formation of a 
glueball sooner than the formation of a qq meson. On 
the basis of this observation, we propose a notion of 
gluon-glueball duality. We apply this to comment on cur- 
rent experimental searches for glueballs and to suggest a 
method that could be useful for these searches. An out- 
line of the paper is as follows. In Sect. II we review 
quark-hadron duality. In Sects. Ill and IV we give some 
background on glueball properties and searches. In Sect. 



V we introduce the notion of gluon-glueball duality in 
two simplified contexts, and in Sect. VI we discuss it in 
full QCD. We point out that in studying the production 
and decay of glueballs, it is useful to analyze the tempo- 
ral evolution of the gg subsystem as it is produced, binds 
to form a proto-glucball, mixes with qq components, and 
finally decays. Section VII suggests some future lattice 
gauge measurements that are relevant to gluon-glueball 
duality, while in Sect. VIII we apply our observations to 
experimental searches for glueballs. 



II. QUARK-HADRON DUALITY 

We first give some background on ideas of duality in 
hadronic physics. The reader who is familiar with this 
material can skip this section and proceed directly to 
our new observations in Sects. V and VI. The idea of 
quark-hadron duality in several related forms @, 0] dates 
back to the early period in the development of the quark- 
parton model. In the Bloom-Gilman form Q, it states, 
roughly speaking, that in a reaction such as a electron 
scattering off a nucleon, the sum of the cross sections for 
the full set of exclusive hadronic final states Xh that are 
kinematically accessible at a given center-of-mass energy 
E cm = y/s is equivalent to the cross section for the ele- 
mentary reaction e + q — > e + q involving the quarks in the 
nucleon. A similar duality relation applies to charged- 
current neutrino reactions such as + N — ► fi + X^. Let 
us denote the four-momenta of the incident and scat- 
tered leptons as l\ and £2, with ii — £2 = q, (^i)iafe = E, 
{i^iab — E' , and the four-momenta of the target nu- 
cleon and final hadronic state as p and px ■ We further 
recall the standard Bjorkcn variables x — ~ q 2 /(2q ■ p) 
and y = q ■ p/t\ ■ p = (E — E')/E. Then this duality is 
the statement that 
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where / denotes all of the charge —1/3 quarks (and 
charge —2/3 antiquarks) that can participate in this re- 
action. At a fundamental level, this duality is justified by 
the asymptotic freedom of QCD Q . In the deep inelastic 
scattering of an electron or neutrino off of a nucleon N, 
the hadronic part of the cross section involves the tensor 

W^(q,p) = ^(7V|J, I |X)(X|jt|iV)(27r) 3 ^ + g-px) 

x 

d 4 z 



« / -^e~^(N\J,(z)M0y\N) , 



(2.2) 



where X denotes a hadronic final state and J ^ is the 
respective electromagnetic or weak (charged or neutral) 
current. One then uses the Wilson operator product 
expansion to express the bilocal product of currents in 
terms of a sum of local operators, applicable near to 
the light cone z 2 — > 0, as enforced by the kinematic 
conditions — q 2 >> Xq CD and q ■ p >> Xq CD , where 
Aqcd — 300 MeV is the scale where QCD confines and 
spontaneously breaks chiral symmetry. This enables one 
to express the deep inelastic scattering off the nucleon 
in terms of the scattering off of quarks. The asymptotic 
freedom of QCD has the consequence that these quarks 
are quasi- free when probed at short distances. Similarly, 
away from particle thresholds, one can calculate the to- 
tal cross section for e + e~ — > hadrons at center-of-mass 
energy y/s in terms of the cross section for e + e~ — > qq, 
% y/s. One can consider the cross section for 
hadrons, smeared over resonances, to be equiva- 

summed over 



where 2m q < ^/s 



lent to the elementary reaction e + e —> QQ, 
the kinematically accessible quarks [9|-(llj 
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In the full QCD theory, the notion of quark-hadron 
duality is naturally generalized to parton-hadron dual- 
ity, where the partons include both quarks and gluons, 
and the hadrons are understood to include not only qqq 
baryons and qq mesons, but also hadronic mass eigen- 
states that are linear combinations of qq and gg, ggg, 
etc. Possible exotic color-singlet hadrons such as, in the 
bosonic sector, qqqq and qqg can, in principle, also be 
included in this set of physical states. In one sense, 
this duality amounts to the statements that (i) there 
is a complete orthonormal basis of perturbative quark 
and gluon states forming the Fock space of perturbative 
QCD, and there is a complete orthonormal basis of phys- 
ical color-singlet hadronic mass eigenstates forming an- 
other Fock space; and (ii), given the asymptotic freedom 
of QCD, the cross section for an inclusive reaction involv- 
ing the contributions of many exclusive physical channels 
with smearing over resonances as appropriate, can be ex- 
pressed in terms of the corresponding cross section in 
terms of the elementary partonic degrees of freedom. In 
another sense, one can think of it as somewhat analo- 
gous to a Mittag-Leffler expansion, in which a function 



is written as a sum over its poles. In this context, one 
may recall that the Mittag-Leffler expansion of the Eulcr 
beta function forms part of the mathematical basis of the 
s-t duality in the Veneziano and Virasoro amplitudes in 
hadronic string theory p^|-(l~6j. 

A specific qq <-* meson duality (and the analogous 
gg <->■ glueball duality to be introduced next) is particu- 
larly useful. This is especially the case if one considers 
the large- N c limit of QCD [jj, [H, El- For large N c , 
baryons become very heavy, and the kinematically ac- 
cessible hadronic states Xh directly produced in e + e~ 
annhilation are qq mesons. Since the decay rate of such 
a meson or glueball vanishes in the large- N c limit, meson 
resonances are narrow in this limit. The energy integral 
in Eq. (|2.3|) then becomes essentially a summation over 
the contributions of these resonances. 

Mesons and baryons are observed to lie on approxi- 
mately linear Regge trajectories of the form 



a(m 2 ) = ao + a'm 2 



(2.4) 



with respective intercepts ao and a common Regge slope 
a' = 0.9 GeV -1 . Physical meson states occur where 
the angular momentum a (to 2 ) is equal to a non- negative 
integer. This behavior was originally motivated by anal- 
ysis of potential scattering and was elegantly explained 
by hadronic string theory (the dual resonance model), 
according to which a meson is a mass eigenstate of an 
open string. It is believed (although it has not been 
proved) that the large-iV c limit of SU{N C ) QCD re- 
produces features of a hadronic string theory. In the 
hadronic string model, the string tension a = l/(27ra'), 
so that y/a ~ 0.42 GeV. Physically, this string tension 
represents the energy per unit length of the chromo- 
electric flux tube between the q and q forming the me- 
son. An example of a Regge trajectory is that for the 
5=1,7=1 (isovector) mesons, which includes p(770), 
a2(1320), p3(1690), and ci4(2040), with increasing values 
of J indicated as subscripts (where J = L + S). The 
radial excitations p' — p(1450), p" = /?(1700), etc. are 
on so-called daughter trajectories, forming a horizontal 
line in the plane with horizontal and vertical axes corre- 
sponding to s = to 2 and J, respectively. 

A feature predicted by hadronic string theory (predat- 
ing QCD) and consistent with data is that the density of 
qq meson states as a function of mass to grows rapidly 
with to. This is also the case for a specific flavor state 
such as ud and specific values of J, parity and charge 
conjugation quantum numbers, such as J PC = 1 . Let 
us denote the density of meson states, i.e., the number of 
states at a given mass to, counting those on the leading 
and daughter meson trajectories, as 



n(m) 



A I 



dn M {m) 
dm 



(2.5) 



where M stands for "meson". For the (bosonic) string 
in d spacetime dimensions, the meson density of states 
n{m)M, grows exponentially fast for to 2 >> (a') -1 (20l . 
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III. REMARKS ON GLUEBALL PROPERTIES 



M 



-(d+l)/4 



exp 



irmyj(2/3)(d-2)a' (2.6) 



where d is the spacetime dimension. Hence, at sufficiently 
high mass, these resonances overlap. Indeed, even before 
one takes account of this asymptotic exponential growth 
in the density of states, the hadronic string model already 
implies that they will overlap, because on the leading 
Regge trajectory, Eq. (|2.4p shows that two successive 
meson states with the same J , that differ by two units 
of L and J, satisfy AJ = 2 = a'(m 2 L+2 — m\), so that 



m L+2 ~ m L 



a'(m L +2 + m L ) 



(2.7) 



Hence, as the masses of these states increase, their mass 
difference decreases, and eventually becomes less than 
their widths, so that they overlap. This happens when 
the mass difference to_l+2 ~ Wj, becomes comparable to 
to either of the widths or Tl + 2- For the present illus- 
trative purposes, we approximate these as being roughly 
equal, and denote them as T, which we take to be 
r ~ 0.25 GeV. Setting to_l+2 — tul = T and solving, 
we get 



m L +2 + niL 



1 



4.5 GeV 



(2.8) 



Thus, as (light-quark) meson masses increase beyond this 
scale, the states in their spectrum tend to merge. In the 
upper end of the mass region of interest here, from about 
1.5 to 3 GeV, the asymptotic condition m 2 >> (a') _1 
begins to be satisfied, so the formula (|2.6p is relevant. 
A hadronization model based on the chromoelectric flux 
tube between a q and q in conjunction with a Schwinger 
mechanism was given in Ref. |22j. The non-Abelian 
generalization, in which a constant chromoelectric field 
creates gluo ns, was analyzed in terms of relevant invari- 
ants in [23l . |24| . The flux-tube mechanism is incorpo- 
rated in current hadronization computer programs such 
as PYTHIA (25[. Because of the increasing density of 
meson states for masses m 2 >> (a') _1 , the cross section 
for e~e + — > qq — > hadrons then becomes a continous 
curve which, according to the duality assumption, coin- 
cides with the continous perturbative curve. 

An important feature concerns the behavior in the 
mass region below approximately 3 GeV. The asymptotic 
freedom and precocious scaling properties of QCD make 
quark-hadron duality a property that is effectively local 
in mass already at masses that are only modestly greater 
than Aqcd- Thus, the p and p' of masses 0.77 and 1.45 
GeV largely account for the contributions in their mass 
region to finite-energy sum rules [12|. This is also man- 
ifest in Bloom-Gilman duality [7j. A difference is that 
Dolen-Horn-Schmid duality applies to 2 — > 2 reactions 
involving onshell hadrons, e.g., 7r + 7r~ — > 7r + 7r~. Simi- 
larly, Bloom-Gilman duality applies to reactions such as 
exclusive electroproduction, e.g., e + p — > e + p + ir°. 



In this section we note some properties of glueballs that 
we will use in our analysis. An especially important and 
relevant property that motivates our new suggestion is 
the density of states, but we begin with some basic facts. 
Since the gluons are bosons, Bose statistics implies that 
the total glueball wavefunction is symmetric under in- 
terchange of any two gluons. A difference between qq 
mesons and glueballs is that although a confined quark 
picks up a (gauge-invariant) dynamical, constituent mass 
of order Aqcd, one cannot ascribe a mass in the same 
manner to a bound gluon, since this would violate the 
color gauge invariance. This means that while a con- 
stituent quark model can provide a good description of 
baryons and qq mesons (see, e.g., [26| for a recent dis- 
cussion and references to the literature), one cannot de- 
scribe the glueball in quite so simple a manner. Further- 
more, in the time evolution of an initial gluonic state, 
the splittings g — > gg can occur in a manner that is lead- 
ing in 1/N C , in the large- N c limit. This is different from 
the time evolution of a qq state, for which the transition 
q — > q + g is suppressed in the large- N c limit. Thus, here 
a physical state denoted as gg strictly refers only to a 
state whose quantum numbers are most simply attainable 
via a (color-singlet) combination of two gluons. Keep- 
ing this caveat in mind, the lowest-lying glueballs can be 
modelled as gg bound states. For these, in the Clebsch- 
Gordon decomposition of the gg SU(3) C representations 
8x8, the singlet appears as a symmetric combination. 
Hence, the product of the space and spin wavefunctions 
must be even under this interchange. The spin wavefunc- 
tion involves the addition of two spin-1 angular momenta. 
If the resultant spin of the gg combination is S — or 
5 = 2, this spin wavefunction is even, so the relative 
angular momentum must also be even, and the ground 
state is L = 0. With P = {-l) L and C = (-I) i+S for 
this combination of two bosons, one thus expects that 
the lowest two glueball states have (i) S = L = J = 0, 
whence J PC = ++ and (ii) L = 0, S = J = 2, whence 
J PC = 2 ++ . The higher- lying glueball states can involve 
both nonzero internal angular momenta and radial exci- 
tations. 

Estimates of glueball masses and widths have been 
made on the basis of a number of different methods 
PH2]> [IjQj (27|- f57Tl . Continuum approaches include the 
MIT bag model [HI]- [H, flux-tube models, AdS/CFT 
approaches, and calculations based on the Bethe-Salpeter 
equation. Lattice calculations have achieved a rather 
high level of precision (2?j], [4!|-[57|]. These naturally give 
the mass of a particular glueball in terms of the square 
root of the string tension, y/a = 0.42 GeV. For masses 
of glueballs in purely gluonic QCD, recent lattice calcu- 
lations [El|, [HHi 1H3 yield 



m(0 ++ ) ~ 1.7 GeV 



m(2 ++ ) ~ 2.4 GeV 



(3.1) 



(3.2) 
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and 



m((T+) ~ m(0 ++ ') ~ 2.6 GeV , (3.3) 



m(l+ - ) ~ m(2 -+ ) ~ 3.0 GeV , (3.4) 



up to approximately 3 GeV. Here, the ++ glueball is 
a radial excitation of the ++ glueball. Lattice measure- 
ments of higher-lying glueball masses have been made up 
to roughly 5 GeV [5lll55l. . Some unquenched calcu- 
lations have also been reported [56| . 

In the context of the Regge or hadronic string model, 
glueballs correspond to closed strings, which have a 
Regge slope equal to half of the Regge slope for open 
strings: 



The relations (|3 - 9[) and (|3 . 10[) follow from direct diagram- 
matic 1/N C counting. As expected from the close corre- 
spondence between the \aige-N c limit of QCD and the 
hadronic string picture, they can also easily be under- 
stood in a string picture. The decay of a qq meson res- 
onance (an open string) takes place via a single cut in 
the string (flux-tube), whereas the decay of a glueball re- 
quires a first cut to transform it from the initial closed 
string to an open string and then a second cut to produce 
the two-meson (e.g., nn) final state. With each cut being 
suppressed by a 1/N C factor, the results on Tm and T G b 
follow. Reverting from the large- N c limit to real QCD, 
actual estimates of glueball widths have varied widely, 
ranging from a few MeV to O(10 2 ) MeV 1^3, iHE El • 



(a') G B = y • 



(3.5) 



It follows that for m 2 >> (a') -1 , the density of states 
for glueballs (closed strings), n(m) G B = dn(rn)GB /dm, 
is exponentially smaller than the density of states for qq 
mesons (open strings), u(to)m- Quantitatively, from Eq. 
PH]) and Eq. (|3~5)) . one finds that, for m 2 » (a')' 1 , 
the ratio of these densities of states is 

n(m) GB _ 9 -(d+i)/4 



exp 



n{m)u 

With a' = 0.9 GeV- 1 (and d 
n(m) GB 



(v^-l) A /(d-2)a'/3 

(3.6) 



n(m) 



4), 

0.3 for m 



M 



and 



n(m) 



CB 



n(m) 



0.1 for 



2 GeV 



3 GeV 



M 



(3.7) 



(3.8) 



Thus, Eq. (|3.6[) indicates that the spectrum of low-lying 
glueball states is more sparse than that of the isoscalar qq 
mesons in the mass region from 1.5 to 3 GeV. To within 
the theoretical and experimental uncertainties, this is 
consistent with the data: for example, in the mass region 
1.3 to 2 GeV, there are the following scalar ++ states 
/o(1370), /o(1500), /o(1710), and indications from recent 
BES data of an / (1790) and / (1810) Q. The lattice 
estimates (to be discussed next) indicate that in this in- 
terval of masses, there is one ++ glueball expected. In 
this channel, this gives a ratio of n{m) G B /n{m)M ~ 0.25. 

Estimates have also been made of glueball widths. In 
the limit N c — > oo with g 2 N c fixed and finite [13, EH, 
where g s is the SU(3) C gauge coupling, the width of a 
glueball vanishes like 



r 



1 



A 



GB 



QCD 



T GB N 2 

while the width of a qq meson M vanishes like 



1 



A 



M 



OC 



QCD 



(3.9) 



(3.10) 



IV. PREVIOUS SEARCHES FOR GLUEBALLS 

Here we briefly review results of previous searches for 
glueballs. There is an extended literature dealing with 
search criteria and analysis of data [E Q , 0, [j| , [H 



87]. 



One signature is that glueballs would not fit into the stan- 
dard set of qq states, including their angular momentum 
and radial excitations. Second, since the gluons carry no 
electric charge, one expects a small branching ratio of 
glueballs into photons. Third, since the gluons carry no 
flavor, it was originally expected that the decays of these 
states should be flavor- independent, up to phase space 
considerations. On the other hand, however, it has been 
suggested that for J — glueballs, there should be helic- 
ity suppression of decays to light-quark hadrons, at least 
if the decay amplitude element can be accurately mod- 
elled beginning with emission of a single qq pair [6l| ; if it 
involves higher initial multiplicity of (anti)quarks, then 
this helicity suppression would be reduced [621 ] . Fourth, 
some glueball states have exotic values of J PC that can- 
not be obtained from qq. 

Experimental searches for glueballs have been car- 
ried out at many laboratories. Experiments using e + e~ 
annihilation include Mark III and the Crystal Ball at 
SPEAR, the subsequent Crystal Ball experiment at 
DORIS, and experiments at other laboratories, including 
Orsay CESR, Novosibirsk, BES, BABAR, and Belle [J- 
11, M~ |Z3- We focus first ° n the isoscalar, J PC = 0++ 
channel, since the lightest pure glueball has these quan- 
tum numbers. There are three prominent isoscalar, 
Lorentz scalar ++ meson resonances between about 
1.0 and 1.7 GeV, namely the /o(1370), / (1500), and 
/o(1710). The quark model is only expected to produce 
two such states, which would have S — 1, L = 1, J = 
and be the analogues of the flavor SU(3) octet and singlet 
pseudoscalar mesons, r/ and rj . The fact that there are 
three fo states in this range is thus one of several pieces 
of evidence suggesting that the third may be primarily 
a glueball. The / (1370) is quite broad, with T - 300 
MeV, while the / (1500) and / (1710) have widths of 
roughly 100-140 MeV More recently, The Beijing 
e + e~ collider BES has found evidence for an /q(1790) 
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and /o(1810) [![. Several theoretical fits to these data 
have been performed 0, Q • The authors of some of these 
fits concluded that the lightest glueball forms a primary 
component in the / (1500) [zl, M HH HI , while others 
concluded that this lightest glueball forms the primary 
component in the /o(1710) [5l|, [83| and still others in- 
voked important contributions from qqqq states |85| 
(see also [34|). Further data and analyses should help to 
elucidate this situation [l|, 0, [86| . 

V. GLUON-GLUEBALL DUALITY IN TWO 
SIMPLIFIED CONTEXTS 

To explain our notion of gluon-glueball duality, we be- 
gin with two simplified forms of QCD, namely (i) with- 
hout any light quarks, and (ii) in the large-iV c limit. Let 
us first consider the case of no light quarks. For defi- 
niteness, we imagine the standard model with one gen- 
eration of fermions with quarks U and D having masses 
mu, mD >> Aqcd- We denote these quarks collec- 
tively as Q. We next consider the favored reaction for 
glueball production, namely the production, in e + e~ an- 
nihilation, of the orthoquarkonium QQ state, followed 
by its radiative decay QQ — » jgg. An important fea- 
ture of this world is that a number of the lowest-lying 
glueball states would be stable. Indeed, using the lat- 
tice estimates of low-lying glueball masses listed above, 
all six of the states listed would be stable; in order for 
a heavier glueball to be kinematically allowed to decay 
to two of the lightest glueballs, it would necessarily have 
a mass greater than about 3.4 GeV. Thus, the invariant 
mass distribution dN/dmc for the mass of the gluonic 
states recoiling against the photon in the radiative ortho- 
quarkonium decay QQ — > ^gg, i.e., at the physical level, 
QQ 7 + Xqb, where Xqb denotes a glueball, would 
exhibit very sharp resonances for rrix GB equal to the mass 
of each of the stable glueballs, and then finite-width res- 
onances for the higher-lying unstable glueballs, up to the 
kinematic limit allowed by the mass of the original or- 
thoquarkonium state. The statement of gluon-glueball 
duality would be that, with appropriate smearing, 

/(f) 

J \dmJ GB J \dmj gg 

where the first integral is over physical glueball final 
states and the second integral denotes the perturbative 
calculation of dN/dm, where m is the invariant mass of 
the gg subsystem in the decay QQ — ► jgg. In terms of 
the overall e + e~ cross section, the gluon-glueball duality 
would be the relation, with appropriate smearing, 

da(e+e- -> n 3 Si(QQ) -> 7 + X GB ) 

GB 

~ da(e+e- -> n 3 Si(QQ) jgg) , (5.2) 

where again the second term represents the perturbative 
calculation of the production and decay. 



In the N c — > 00 limit, qq mesons and glueballs be- 
come stable, as indicated by Eqs. (|3.10[) and (|3.9p . Fur- 
thermore, there is no mixing between glueballs and qq 
mesons. Here, gluon-glueball duality takes a particularly 
simple form. With N c large but finite, so as to allow for 
the radiative decay of the heavy orthoquarkonium state, 
this duality would again be expressed via the relations 
(15. ip and (|5.2p . Quark-hadron duality also takes a par- 
ticularly simple form in this large- N c QCD. This type 
of connection between sums over resonances and proper- 
ties of the underlying quarks and gluons was previously 
used with QCD sum rules to study correlators of various 
operators [2§|-[34|. 

VI. GLUON-GLUEBALL DUALITY IN QCD 

We next discuss our notion of gluon-glueball duality in 
real QCD. An important part of our discussion of this du- 
ality in the radiative decay of a heavy orthoquarkonium 
state QQ — » *ygg is a careful treatment of the temporal 
evolution of the gg subsystem, as it is initially produced, 
as the gluons bind to form a proto-glueball, as this glue- 
ball mixes with a qq component, and as it finally decays. 
To understand gluon-glueball duality, it is crucial to an- 
alyze the time evolution and hierarchy of time scales rel- 
evant to the QQ — ► jgg decay, as compared with the 
production of mesons in a reaction such as e + e~ — > qq. 
A general statement concerns the time required for the 
formation of color-singlet states from the respective ini- 
tial qq and gg states. Given the fact that QCD confines 
on a scale Aqcd and that hadrons have a corresponding 
size 

Thad. ~ — ~ 1 fm , (6.1) 

and given the causality condition that information can- 
not be communicated any faster than at the speed of 
light, it follows that a minimum time associated with the 
formation of color-singlet hadronic states is 

thad. = — 0.3 x 1CT 23 sec (6.2) 

c 

(where we have explicitly indicated the speed of light, c). 
This is a rough estimate, accurate to a factor of order 
unity. For example, given that a glueball is represented 
by a closed string, one could consider a special case in 
which the closed string forms a circle, and one might 
argue that it is the circumference of this circle rather than 
the radius that is of order 1 fm. In this case, the radius 
would be l/(27r) fm and the time taken for the formation, 
involving motion of the gluons outward from the center 
of the circle would be smaller than the value given in 
Eq. (|6.2[) by the factor 2ir. Because of the asymptotic 
freedom of QCD, for both (i) e + e~ — * qq at center-of- 
mass energy \fs >> A^qcd and (ii) the radiative decay 
of heavy orthoquarkonium QQ — > ^gg, there exists a 
sufficiently short time t pert . such that for times t < t per t. 
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the physics can be described using perturbative QCD. 
This satisfies the inequality 



has, roughly, 



Lpert. 



< 



1 



A 



QCD 



thad. • 



(6.3) 



Given the precocious scaling behavior of QCD, it is not 
necessary that t pert . << thud.- For the two specifc 
cases under discussion, one could take t 
for the reaction e + e" 
for the decay QQ — > 



pert. 



1 



pert. 



— > qq and t pert , — 1/(2toq) 
jgg. A typical value would be 
1/(3 GeV) ~ 10 -25 sec. Thus, in the reaction 
e ■ e — > gg, during the time interval < i < t per t., 
the g and g recede from each other in an approximately 
perturbative manner, with the first modification being 
the emission of a gluon, leading to a qq subsystem in a 
color octet state together with the emitted gluon g. In 
the radiative decay of the heavy orthoquarkonium state 
QQ ~ * 7ff<7> during the time interval < t < t pert . , the gg 
final-state subsystem mainly evolves into more gluons via 
g — > gg splittings. As noted above, this gluon splitting 
occurs at leading order in the large- N c limit, in contrast 
to the q —* q + g or g — > (qq)s processes, which start to 
mix gg with the initially purely gluonic gg subsystem. 

After a time tMF, where MF stands for "meson for- 
mation" , the initial gg system will bind to form a meson, 
and after a corresponding time taBF, where GBF stands 
for "glueball formation" , the initial gg system will bind 
to form a glueball. From the causality argument above, 
one has the general inequalities 



tMF , tGBF > thad. 



(6.4) 



and hence also the obvious inequalities tMF, tc b f > 
t pert .. The q and q in the meson, and the gluons in the 
glueball have minimum bound-state momenta k m in ~ 
Aqcd because of confinement [35j. Several factors are 
relevant for the hadronic formation times tMF and tcBF, 
including (i) the intrinsic QCD hadronization time scale 
thad. , (h) the mixing of gg and gluonic states to form mass 
eigenstates, (iii) the decay widths I\ of various mesons 
and glueballs and (iv) especially importantly for our cur- 
rent discussion, the density of meson and glueball states, 
n(m)M and n(m)aB- The quantum-mechanical uncer- 
tainty relation AE At > fi/2 implies that the observa- 
tion time interval At needed for an observer to measure 
the spectrum of states with a resolution in mass Am is 
bounded below by At > (fi/2)/Am. Here Am is set by a 
combination of the density of states with the same quan- 
tum numbers (isospin and J ) and by the widths of 
these states. Let us consider a glueball search conducted 
in the range of masses togs = 1-5 — 3 GeV. Given the in- 
equality in the density of glueball versus gg meson states 
in Eq. (|3.6[) . it follows that the time needed to exper- 
imentally measure and resolve glueball states is shorter 
than that needed for gg meson states. Using the hadronic 
string model as a theoretical guide, which is consistent 
with the observed states in the relevant mass region, one 



t G BF 



tMF 



(6.5) 



This leads us to suggest a different picture of glueball 
production than the one that is often used in analyses 
of experimental data on glueball searches. Conventional 
analyses use meson mass eigenstates that are linear com- 
binations of gg states and gluonic states. Our new point 
is that it is crucial to take into account the actual tempo- 
ral formation of the glueball states. Given that the glue- 
ball formation time is shorter than the meson formation 
time, with tcBF — tMF / 4 being a reasonable estimate, 
the glueball forms before significant mixing with the gg 
sector takes place. A concrete realization of both the 
gg <-> meson and gg <-» glueball dualities can be obtained 
as follows. Starting with an initial pure gg entrance state, 
we implement the duality by letting the unitary QCD 
evolution operator, formally expressed as U(t) = e~ tHt , 
operate on this state, where here H denotes the QCD 
Hamiltonian, yielding 



U{t M F) \qq(t = 0)) = \M, meson) 



(6.6) 



That is, the evolution over this time interval will yield a 
physical gg meson resonance. In a similar manner, in a 
purely gluonic sector 



U{t GBF ) \gg(t = 0)) = \GB, glueball) 



(6.7) 



The two gluons in the gg subsystem produced in the ra- 
diative orthoquarkonium decay QQ — > jgg emerge from 
spacetime points that are separated by a small distance 
Ar ~ 1/ttiq, where Q = c or b is a heavy quark. This 
is not precisely the same as the production of a scalar 
glueball by the action of the local operator 



S{x)=G liV {x)G^ v {x) 



(6.8) 



on the vacuum. However, a semiclassical argument leads 
to the conclusion that the gg usually bind with L = 
relative orbital angular momentum. For example, in the 
case Q = b, the spatial separation of the points where 
the two gluons are emitted is Ar ~ 1/mt ~ 0.2 GeV" 1 . 
The 3-momenta of the gluons in the gg rest frame are 
\k g \ ~ mcs/2 ~ 0(1) GeV. The resultant average value 
of the relative orbital angular momentum is 



(L) 



\k g \Ar < 0.2 



(6.9) 



Hence, one expects that this production mechanism will 
yield mainly glueball states with L = 0, namely the ++ 
and 2 ++ mentioned before. 

In the radiative QQ — * jgg decay of orthoquarko- 
nium, the gg subsystem is manifestly purely gluonic to 
start with, and mixing with gg components occurs sub- 
sequently. In the large- N c limit, this mixing is sup- 
pressed by 1/N C , which has led to the common expec- 
tation that there could be hadrons that are primarily 
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gluonic, with only a small qq component. Our estimate 
that tcBF — ^mI^i m conjunction with the suggestion 
from large- N c arguments that mixing of gluonic and qq 
components may be rather small, leads us to the impor- 
tant inference that the gluon-glueball duality could hold 
reasonably well in full QCD as well as in the simplified 
contexts which we initially used to introduce it. It is un- 
derstood that there will be some corrections due to the 
the mixing of the gluons in a primarily gluonic hadron 
with qq states. 

Let us next consider the longer times required for the 
qq meson, or the glueball, to decay into hadrons that 
are stable with respect to the strong interactions. The 
formation and decay times for the meson resonances are 
comparable, although tm <^ tMF > t pe rt.> an d similarly 
for the glueballs, one has tqb ^ ^gbf > t pert .. This 
is to be contrasted with the situation for a very heavy 
quark, namely the top quark, which decays weakly be- 
fore it can form color-singlet hadronic tt or tq states. To 
the extent that the large- N c limit is applicable to QCD, 
one expects that, other factors such as phase space being 
equal, the lifetime tqb might be somewhat longer than 
tm, i.e., the glueball width might be somewhat smaller 
than that for a qq meson of comparable mass. However, 
in actual QCD, glueball widths may not be suppressed, 
and may, indeed, be of order 100-300 MeV. This would 
be somewhat analogous to the situation with the 77' me- 
son; in the iV c — + 00 limit (with A = g N c fixed), in- 
stanton effects are exponentially suppressed by the fac- 
tor exp(— 8n 2 /g 2 ) = exp(— 87T 2 iV c /A), so that \J(1)a is a 
good global symmetry and the isoscalar pseudoscalar me- 
son r{ is an approximate Nambu-Goldstone boson. How- 
ever, in real QCD the 7/ is rather heavy, with a mass 
of 958 MeV. An important point is that, with the hi- 
erarchy of time scales that we have noted, the glueball 
decays by popping two pairs of light qq quarks out of the 
vacuum to produce the two final-state mesons (7r's, K's, 
etc.). This process is essentially equivalent to the process 
by which the initially purely gluonic state acquires a qq 
component. 



VII. FURTHER POSSIBLE INSIGHT FROM 
LATTICE QCD 

Lattice calculations have the appeal of providing a 
fully nonperturbative tool for studying the properties of 
QCD, and the advantage of being able to be continu- 
ally improved with the use of larger lattices, longer run- 
ning times, improved lattice actions, and careful analysis 
of statistical and systematic uncertainties. Most lattice 
QCD calculations of glueball masses have been performed 
using the quenched approximation. Some unqucnchcd 
calculations have also been reported (56|. Both the ne- 
cessity of evaluating the fermion determinant and the 
related presence of disconnected flavor loops appearing 
in unquenched calculations make these calculations more 
difficult than computations in quenched QCD. We sug- 



gest that it would be worthwhile for lattice gauge simu- 
lations to address some of the issues that we have raised 
in this paper. We are interested not just in minor shifts 
of the glueball spectrum, but rather in finding the time 
t* by which the admixture of the initial glue state with 
the qq and multi-quark states becomes significant. For 
this purpose it could be useful to study the correlator 
C(t) — (S(0)S(t)) of the above-mentioned scalar glue- 
ball operator and examine how its Euclidean time de- 
pendence might differ from a simple exponential of the 
form exp(— m(0 ++ ) t). (Here, it is understood that one 
would ideally have removed the effects of higher-lying 
glueball states with the same J PC = ++ and also that 
one would have taken account of effects due to periodic 
lattice boundary conditions.) For long, asymptotic times 
t such that t >> l/(2m jr ), the behavior of this scalar 
correlator C(t) is controlled by the lowest-mass s-channel 
threshold, namely that for the 2ir final state, but we are 
interested in shorter times. Similar calculations could be 
performed for the 2 ++ glueball state by using an appro- 
priate color-singlet tensor correlator. Assessing the full 
lifetime until the glueball decays into final hadrons that 
are stable with respect to the strong interactions is chal- 
lenging, but is not essential for our purposes here. 



VIII. APPLICATION TO EXPERIMENTAL 
SEARCHES FOR GLUEBALLS 

In this section we apply our notion of gluon-glueball 
duality to suggest a method that could be useful in exper- 
imental searches for glueballs in radiative orthoquarko- 
nium decays, in particular, those involving the T(15), 
T(2S), and T(3S) states. There are very high-statistics 
data sets from radiative J/ip decays, which have been 
used quite effectively for glueball searches. However, ra- 
diative T decays allow one to search in a wider mass 
range and reduce phase space suppression for decays into 
final states involving more massive glueballs. While a 
major purpose of the experiments at BABAR and Belle 
was to study B physics and CP violation, they accumu- 
lated of order 10 9 events from decays of T(1S), T(25), 
and T(35), as well as the T(4S) state that provided a 
copious source of Bd mesons [70, [z3- These data ex- 
tended the already impressive data sets collected by the 
CLEO experiment at CESR in its later years of high- 
intensity running [r35l . Ir36j . For an average radiative de- 
cay branching ratio of 1.5 % we expect of order 1.5 x 10 7 
radiative decays in the BABAR and Belle data. In the 
radiative T(nS) — > jgg (with n — 1,2,3) we label the 
four-momenta of the outgoing photon and gluons as fc 7 , 
k±, and &2, and recall that it is necessary to symmetrize 
the amplitude under the interchange k\ <-> k% to take ac- 
count of the two identical bosons in the (perturbative) 
final state. If one makes the approximation, in the per- 
turbative calculation of the amplitude, that the outgoing 
(massless) gluons interact only very weakly with each 
other, it follows that the three invariant mass combina- 
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tions (fci + fc2) 2 , (fc 7 + fci) 2 , and (fc 7 + fc^) 2 are uniformly 
distributed over the Dalitz plot, which becomes an equi- 
lateral triangle. In this Dalitz plot, the region of interest, 
which is assumed here to be dominanted by the lowest- 
lying glueballs, is then a rectangular strip adjacent to the 
bottom of the triangle. The total area of this region is 
2 • (2.7) 2 /100, i.e. 15 % of the total area of the Dalitz 
plot and hence includes approximately 2 x 10 6 events. 
The notion of gluon-glueball duality that we have dis- 
cussed then leads us to the suggestion to analyze the 
inclusive mass distribution of these 2 x 10 6 events. This 
avoids any bias due to post-selection by the final chan- 
nel (which might prefer specific final qq or multiquark 
resonances). Our use of gluon-glueball duality is analo- 
gous to the use of quark-hadron duality in the sense that 
both of these dualities relate inclusive channels and sums 
of exclusive channels in the respective particle processes. 
For notational simplicity, we denote X = Xqb- Clearly, 
only a crude resolution AMj ~ 0.5 GeV is needed to re- 
solve the two well-separated lowest-lying glueball states 
with J PC = 0++ and 2++ (or the excited 0++ state). 
Let us denote the invariant mass squared of the gg sub- 
system as M x — (fci + /C2) 2 and take particle energies to 
be measured in the rest frame of the decaying QQ ortho- 
quarkonium state. The elementary kinematic relation 

M QQ = ( fc 7 + ki + k 2 ) 2 = 2E y {M QQ - Ey) + M x (8.1) 



implies that 



AM X = 



(2£ 7 - Mqq)AE~i 
M x 



(8.2) 



As an illustration, we consider the BABAR detector [70j; 
similar numbers apply for the Belle detector [7lj . The 
fractional resolution (AE 7 )/ E 7 , of the measurement of 
the photon energy by the electromagnetic calorimeter of 
this detector varies from from about 2 to 3 % over the 
range of £ 7 from - 8 GeV to 1 GeV Hence, the 

resultant resolution AMx from Eq. (|8.2|) , for the ra- 
diative decay of the T(IS'), varies from approximately 
0.55 GeV to 0.27 GeV as M x varies from 1.7 GeV to 2.4 
GeV. For the radiative decay of the T(25) the resolution 
AM X varies from about 0.65 GeV to 0.33 GeV as M x 
varies from 1.7 GeV to 2.4 GeV. Considering the very 
high statistics of the data sets obtained by BABAR and 
Belle, this analysis could give useful information about 
glueballs via broad deviations from the phase space dis- 
tribution that would occur in their absence. This analysis 



presumes that one takes careful account of pure quantum 
electrodynamic (QED) backgrounds and corrections. By 
insisting on some hadronic activity in the detector, one 
may reduce such QED backgrounds without excessive bi- 
asing such as would result if one were to fully reconstruct 
the final hadronic state. Obviously, the experimental 
procedure sketched here in the broadest terms is chal- 
lenging. Nevertheless, one has observed how much useful 
new data BABAR and Belle have obtained concerning 
new hadronic states involving charm quarks, including 
A(3782), new D s states, and others. Provided that our 
analysis of the time evolution of the glueball production 
process discssed above is correct, then we believe that 
these facilities have the potential to considerably clarify 
the lingering puzzles in glueball physics. 



IX. CONCLUSIONS 

In this paper we have presented a different picture of 
glueball production than the one commonly used in cur- 
rent analyses of data. Using the hadronic string model, 
we have given a quantitative estimate of the smaller 
density of states of glueballs (closed strings) in the re- 
gion of ~ 2 GeV, as compared with qq mesons (open 
strings), and, from basic quantum mechanics, we have in- 
ferred a resultant hierarchy of formation times of observ- 
able (resolvable) glueballs, as compared with qq mesons, 
namely Eq. (|6.5p . On the basis of this, together with 
the suggestion from the large--/V c expansion that mix- 
ing between glueball and qq states may be suppressed, 
we have argued that the glueballs produced in radia- 
tive orthoquarkonium decay could plausibly form with- 
out substantial mixing with qq states. This motivates a 
notion of gluon-glueball duality, which we have presented, 
namely that the summation over sufficiently many glue- 
ball states produced in radiative orthoquarkonium decay 
QQ ^ 199i appropriately smeared, could be well fit with 
the perturbative calculation of this process. We have ap- 
plied this notion of gluon-glueball duality to suggest a 
method that could be useful in experimental searches for 
glueballs using radiative decays of the T(15), T(25), and 
T(35) states using the large data sets that are currently 
available on these decays. 
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